DDC  _fjlej)opyj  WD  A 0 6 9 8 8 4 


AFOSR-TR-  7 9 - 0 67  1 


ENGINEERING 

DESIGN 

CENTER 


SCHOOL  OF 
ENGINEERING 

NORTH  CAROLINA 
STATE  UNIVERSITY 

RALEIGH 

NORTH  CAROLINA 


E?  PTTV' ,r’  rl 

H f.  y t 

& !■  Iff  I 


EDC-79-1 


An  Actuator  Disk  Analysis  of  an 
Isolated  Rotor  with  Distorted  Inflow 


DDC 


J.  N.  Perkins 


ifscmiiEjnj 

JUN  13  1979  | 


r liULb^ 


Report  on  AFOSR  Contract  F44620-76-C-0055 
for  a cooperative  program  between 

United  Technologies  Research  Center 

and 

North  Carolina  State  University 


releag®? 

Approvod  for  ^UBited# 

distribution 


I 


EDC-79-1 


i 


An  Actuator  Disk  Analysis  of  an 
Isolated  Rotor  with  Distorted  Inflow 


by 

J.  N.  Perkins 

Report  on  AFOSR  Contract  F44620-76-C-0055  for  a 
cooperative  program 

between 


United  Technologies  Research  Center 
East  Hartford,  Conn. 

and 


North  Carolina  State  University 
Raleigh,  North  Carolina 


Engineering  Design  Center 
School  of  Engineering 
North  Carolina  State  University 
Raleigh,  N.  C.  27650 

March  31,  1979 


AIR  FORCE  OFFICE  0?  SCIENTIFIC  RESEARCH  (AFSC) 

NOTICE  OF  TRANSMITTAL  TO  PPG 
This  technical  i .r-  . t i rtfvie*ed  and  is 

approv'd  for  p ■!  1 . lease  1 A>V  AFR  li.'U-12  (7b). 

Distribution  i.»  un  1 ircittu. 


A.  0.  IS  LOSE 

Technical  In/oi  aallon  Officer 


X 


r 


ii 


Abstract 

An  analytical  study  of  the  passage  of  distorted  flow  through  an 
isolated,  high  hub-tip  ratio  axial  compressor  is  reported.  The  analysis 
involves  the  coupling  of  the  unsteady  blade  row  aerodynamic  response  (with 
the  flow  immediately  upstream  of  the  blade  row  prescribed)  with  the  unsteady 
duct  flow  both  upstream  and  downstream  of  the  blade  row  under  the  influence 
of  the  blade  row  loading.  The  numerical  solution  uses  a starting  procedure 
which  does  not  require  the  inlet  plane  to  be  far  enough  upstream  to  be 
unaffected  by  the  presence  of  the  blade  row.  Hence  any  experimentally  de- 
termined distortion  at  any  arbitrary  distance  upstream  of  the  blade  row  can 
be  modelled.  The  results  obtained  indicate  that  the  predicted  pressure 
profile  at  the  blade  row  exit  is  strongly  dependent  on  the  experimentally 
determined  steady-state  loss  and  exit  flow  angle  curves,  but  is  almost 
independent  of  the  magnitude  of  the  first  order  lag  coefficient  used  to 
represent  the  boundary  layer  time  delay  through  the  blade  passages. 
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Introduction 

The  successful  prediction  of  the  performance  of  a bladed  compressor 
requires  adequately  quantitative  mathematical  methods  for  describing 
the  flow  behavior  encountered  in  the  machine  in  terms  of  accurately  measur- 
able quantities.  In  the  present  study,  the  response  of  an  isolated  rotor 
to  inlet  flow  distortions  is  considered. 

Of  the  various  inlet  flow  distortion  modes,  the  uniform  radial 
condition  with  circumferentially  distorted  axial  flow  has  been  found  to  have 
large  effects  on  compressor  performance.  Thus,  this  type  has  been  exten- 
sively investigated  both  analytically  and  experimentally.  A useful  biblio- 
graphy can  be  found  in  Refs.  1 and  2 which  may  serve  as  guides  for  the 
current  state  of  the  art,  and  may  provide  reference  lists  to  the  vast  liter- 
ature already  existing  in  the  field.  Different  models  have  been  developed 
to  predict  the  response  of  compressors  to  circumferential  inlet  distortions 
and  recent  examples  of  the  more  complete  treatments  are  given  in  Refs.  3-5. 

One  of  the  major  difficulties  inherent  in  developing  such  analyses  is 
that  the  rotor  blades,  which  move  through  a spatially  nonuniform  flow  field, 
are  subjected  to  unsteady  flow.  This  unsteady  performance  is  an  essential 
part  of  any  analytical  perdiction  technique  for  distorted  flow.  Although 
some  success  has  been  achieved,  the  models  presently  used  to  describe  this 
unsteady  response  are  still  quite  rudimentary,  and  there  is  considerable 
room  for  improvements.  One  of  the  major  difficulties  here  is  the  lack  of 
detailed  quantitative  knowledge  of  the  unsteady  flow  field  imposed  on  the 
blades.  In  fact,  due  to  this  lack  of  quantitative  experimental  evidence 
for  isolated  rotors,  very  few,  if  any,  direct  comparisons  between  the  existing 
analytical  models  and  experiment  have  been  made.  Therefore,  North  Carolina 
State  University  and  United  Technologies  Research  Center  have  embarked  on  a 
joint  project  under  sponsorship  of  the  Air  Force  Office  of  Scientific  Research 
to  study  analytically  and  experimentally  the  response  of  a compressor  stage 
to  an  inlet  distortion.  The  work  reported  herein  is  the  first  of  several 
planned  analytical  investigations  and  direct  comparisons  to  experiment  of  this 
p hen omen. 
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Semi -actuator  Disk  Analysis 

A nonlinear,  large  disturbance  theory  was  developed  which  couples, 
interactively,  the  flow  through  the  passages  of  an  isolated,  high  hub-tip 
ratio  axial  compressor  blade  row  and  an  axially  distorted  incompressible 
flow  field.  The  blade  row  analysis  is  based  on  the  time-dependent  energy 
equation  developed  in  Ref.  3.  The  flow  field  analysis  involves  the  two- 
dimensional  form  of  the  time-dependent,  incompressible,  viscous  Navier  • 

Stokes  equations.  Coupling  of  the  two  is  accomplished  by  satisfying  the 
boundary  conditions  relating  the  pressure  change  across  the  blade  row  and 
those  requiring  mass  continuity  through  the  blade  row. 

Since  the  flow  is  assumed  to  be  independent  of  radial  position, 
the  analysis  can  be  carried  out  in  the  two-dimensional  cascade  plane.  This 
is  shown  schematically  in  Fig.  1,  which  represents  the  flow  field  ahead  of 
and  behind  an  axial  compressor  blade  row,  unwrapped  tangentially  and  projected 
onto  a two-dimensional  plane.  The  lines  ABC  and  DEF  represent  a common 
axial  trace  in  the  three-dimensional  compressor  system,  and  this  correspon- 
dence imposes  a periododic  boundary  condition  on  the  flow  velocities.  A 
distorted  inlet  flow  is  initialized  along  the  inlet  boundary  AD,  and  a con- 
stant, although  unknown  a priori,  static  exit  pressure  condition  is  imposed 
at  the  exit  boundary  CF.  The  blade  row  represented  by  the  line  BE  is 
modelled  as  a semi -actuator  disk  across  which  there  is  a change  in  the  flow 
direction  as  well  as  a change  in  total  pressure. 

The  flow  field  upstream  and  downstream  of  the  blade  row  is  required 
to  satisfy  the  incompressible  Navier  Stokes  equations  of  motion  without 
assuming  conservation  of  mass  (see  Fig.  1 for  nomenclature), 
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and  the  equation  of  continuity. 


D'  = 


3U 1 . 3V* 


(3) 


3y ' 3y ' 

The- mass  discrepancy  term,  D',  should  ideally  be  zero.  However,  due  to  in- 
accuracies in  the  finite  difference  computations,  this  term,  while  required 
to  be  small,  is  not  zero  and  is  left  in  the  equations  to  eliminate  a non- 
linear instability  in  the  solution. 
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Introducing  the  nondimensional  quantities: 


u = u'/U  , v = v'/U  , p = p'/pU2 


x = x'/a  , y = y'/Jt  , t = t'U A 


eqs.  (1)  - (3)  become 


3u  3U2  _ BUV  _ 3£  + l_,3£u  + 32U  + 1 3D 

3t  3x  " 3y  3x  Re'ax2  3y2  3 3x 

3V  _ _3v£  3UV  3£  + 1_/3ZV  + 32V  + 1 3Dx 

3t  "3y  3x  3y  Re  cix7  ay7  3 3y; 


0 = 1“  + !*  ( 
3x  3y 

By  differentiating  (5)  with  respect  to  x,  (6)  with  respect  to  y,  and 
(7)  with  respect  to  t,  Poisson's  equation  for  the  pressure  can  be  obtained 
as  follows: 


92u  _ 32U2  _ 32UV  _ a£p  . 1 # 33U  + 33U  + 1 32Dn 
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Note  that  for  the  flow  under  consideration,  the  value  of  D should  be 
very  small  before  the  above  equations  can  be  said  to  model  the  flow  accurately. 
The  mass  discrepancy  term  is  retained  for  the  sole  purpose  of  eliminating  a 
nonlinear  instability  in  the  computer  solution. 

By  assuming  that  the  flow  in  the  blade  passage  of  the  cascade  is  like 
channel  flow,  the  unsteady  blade  row  equations  can  be  derived  as  in  Reference  3: 


- P62)  = c sec  ?(-£)r  + |x[uj  + (V+vj)  ] 


(10) 


where  X is  an  unsteady  loss  parameter  defined  as  the  instantaneous  sum  of  the 
gap-averaged  deviation  of  exit  total  pressure  and  the  dissipated  mechanical 
energy  during  an  unsteady  process.  In  a steady-state  process  eq.  (10)  shows 
that  X reduces  to  the  standard  steady-state  loss  coefficient, 
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Since  the  specific  unsteady  behavior  of  X is  not  known,  the  first-order  lag 
model  of  ref.  (3)  is  used: 


x=xss  - 


(12) 


Figure  2 shows  the  stationary  and  blade  reference  frames.  The  time 
derivatives  in  the  blade  reference  frame  are  transferred  to  the  stationary 
frame  by: 
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- Po2)  - c sec  s(f$  - v'0)  + [x[u!'2  + (V*  + Vl*)2] 


and 


y _ v _ T < _ v 'J(X.  ) 

* Ass  ^ at 1 v ay1' 


Since  continuity  across  the  blade  channel  requires 


u*  = u'  = u'  , 
1 2 


then  from  Fig.  2 it  can  be  seen  that 
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cot  p = L 
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= V,+v? 
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Using  Bernoulli's  equation  at  the  inlet  and  exit  plane. 


P<5  = Pi  + \ p[u'2+  (V’+v')2] 
P6  = p'  + \ P[u'2+  (V'+v')2] 

o ? c. 


Introducing  the  additional  dimensionless  quantities 


V = V'/U,  c = c'/fc  . 


- -'U/* 


T = T 


ri  imaiiUftiir 


(14) 


(15) 


(16) 


(17) 


(18) 


Using  eqs.  (17)  and  (18),  eq.  (14)  can  be  written: 


l(p{-p^)=  c'  sec  5 [|p  - V'f£]  + ^J,2[X+(X-l)cot2P1  + cotzp2]  (19) 


(20) 


and  using  eqs.  (4),  the  dimensionless  blade  row  coupling  equations  become: 


Pj  - P2  = c sec  5 [|^  - v|y]  + |ti2[X  + (X-l )cot2Bi  + cot232]  (21) 


where 


X = XSS  - '(lf  - v!f>  <22> 

The  values  of  X$s  and  &2  as  functions  of  and  the  value  of  t must 
be  determined  from  experiment.  Equations  (5)  through  (9),  (21)  and  (22)  are 
the  partial  differential  equations  to  be  solved.  However,  before  a solution 
can  be  attempted,  it  is  necessary  to  have  suitable  boundary  and  initial 
conditions. 

Figure  1 illustrates  the  six  boundaries  of  importance  in  this  problem. 
The  upper  and  lower  boundaries  are  the  easiest  to  handle.  All  conditions  on 
the  lower  boundary  are  equal  to  those  on  the  upper  boundary.  This  is  equi- 
valent to  saying  in  the  case  of  the  actual  rotor  that  conditions  are  periodic 
around  the  circumference. 

It  will  be  assumed  that  the  inlet  is  far  enough  away  from  the  blade 
row  so  that: 


By  applying  continuity,  another  condition  can  be  found: 

9u  _ p 3v 
3x  ay 

Differentiating  (10)  with  respect  to  x yields: 


(23) 


(24) 
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It  will  also  be  assumed  that  the  total  pressure  remains  constant 
with  time  along  the  inlet.  Differentiating  Bernoulli's  equation  with  re- 
spect to  t gives: 

1 / 3U2  . 8 V2 


IE  = _ 1/ 

3t  2^at 


+ if  ) 

at  ' 


(26) 


For  the  outlet, assumptions  (23),  (24),  and  (25)  are  used  again.  These 
conditions  are  applied  in  (5)  to  formulate  the  Neumann  pressure  gradient 
boundary  condition: 


ap  _ au2  _ ajjv  + l_  a2u 
ax  ax  ay  Re  ay7 

3 U 

Note  that  in  the  above  ^ has  been  neglected. 

The  most  complicated  boundary  conditions  occur  at  the  blade  row 
entrance  and  exit.  Assumptions  (23),  (24),  and  (25)  will  be  utilized  at 
both  entrance  and  exit  along  with: 


(27) 


D = 0 , — = - — (28) 

* ax  ay 

The  pressure  gradient  at  entrance  and  exit  can  be  found  as  before: 


I aP\  _ /dUi  _ /dU-\  _ / dUV\  , I /d~U\ 
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The  unsteady  term  must  be  included  here.  From  continuity  across 


r a 2 U \ 

'9^2 


(29) 


(30) 


(31) 


Subtracting  (30)  from  (29)  and  applying  (31)  leads  to  an  expression 
coupling  the  pressure  fields  on  either  side  of  the  blade  row: 


<S>, 


/ 3P\  _ (3U2\  /3U2} 

'3X  2 '3X  2 " '3X  \ 


+ (i^) 


, 3UVx 
l3y  >■ 


(32) 


Initial  conditions  compatible  with  the  above  boundary  conditions  are 
needed  to  start  the  problem.  The  initial  inlet  axial  velocity  and  pressure 
distributions  are  obtained  from  curve  fits  of  experimental  data.  The  initial 
inlet  axial  velocity  distribution  will  be  assumed  constant  in  x up  to  and  in- 
cluding the  blade  row  exit.  The  circumferential  velocity  will  be  zero  in 
the  entire  upstream  region.  From  this  information  values  of  3i  and  corres- 
ponding values  of  e>2  can  be  generated.  This  will  give  the  circumferential 
velocity  distribution  behind  the  blade  row.  From  continuity: 


(— ) = (— ) 
v3x'2  3y  2 


(33) 


As  a starting  point,  it  is  assumed  that  behind  the  blade  row: 


ly.  - P -A(x-x2) 

3X  ' 3X  2 

Again  applying  continuity: 


(34) 


li.  = fly.)  p -a(x-x2) 
3y  ' 3X  2 


(35) 


Also  along  the  center  line  behind  the  blade  row,  it  is  arbitrarily 
assumed  that  initially 

fx  ■ 0 (36> 

Note  that  D equals  zero  everywhere  initially. 

Also: 

~ = 0 (37) 
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The  next  step  in  the  analysis  consists  of  finite  differencing  the 
above  partial  differential  expressions.  The  differencing  technique  used  is 
that  of  Welch,  et.  al . (Ref.  6).  The  results  along  with  the  mesh  setup  are 
given  in  the  appendix. 

The  numerical  procedure  used  to  obtain  the  solution  to  the  finite 
differenced  equations  and  associated  boundary  conditions  as  given  in  the 
appendix  employs  the  algorithm  for  solving  elliptic  problems  by  direct 
marching  methods  given  by  Roach  in  Ref.  (7).  This  method  was  used  because 
it  was  found  to  be  at  least  two  orders  of  magnitude  faster  than  any  other 
method  which  did  not  diverge.  While  the  marching  algorithm  is  field  size 
limited,  no  problems  were  noted  with  the  present  field  size  which  extended 
from  approximately  4 chord  lengths  upstream  of  the  blade  rov  to  approximately 
15  chord  lengths  downstream  of  the  blade  row.  The  starting  procedure  is  to 
initialize  an  axial  velocity  distortion  at  t = 0 at  the  inlet  plane,  AD,  and, 
assuming  no  work  being  done  on  the  fluid  at  this  plane,  maintain  the  resulting 
total  pressure  at  the  inlet  constant  for  all  later  times.  This  does  not  re- 
quire making  the  usual  assumption  that  the  inlet  plane  is  far  enough  up- 
stream to  be  unaffected  by  the  presence  of  the  blade  row  at  BE.  At  each 
time  step  over  t > 0,  the  propagation  through  the  entire  flow  field  of  the 
disturbances  created  by  the  blade  row  is  computed  subject  to  a simultaneous 
satisfaction  of  all  boundary  conditions.  The  computation  is  terminated 
when  a steady  state  solution  is  obtained.  The  resulting  axial  distortion 
is  then  compared  to  the  experimentally  determined  distortion  at  the  up- 
stream location  where  the  measurements  were  made.  If  the  two  distortions  are 
not  identical,  the  initial  distortion  at  the  inlet  plane  is  changed  and  the 
computation  repeated. 

Numerical  Results 

The  finite  difference  equations  and  associated  boundary  conditions 
given  in  the  appendix  were  programmed  for  the  NCSU  IBM  370  digital  computer 
and  a number  of  cases  corresponding  to  the  experimental  results  presented  in 
Ref.  8 were  examined.  The  required  values  of  X$s  and  b2  as  functions  of  $i 
which  were  used  in  the  computations  are  presented  in  Figs.  3 and  4 respectively, 
and  were  obtained  from  the  experimental  results  of  Ref.  8 . As  was  noted  in 
Ref.  8,  both  the  loss  and  turning  characteristics  as  obtained  in  that  study 
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were  thought  to  be  in  error.  Thus,  the  results  obtained  from  the  present 
theoretical  model  are  also  thought  to  be  in  error.  In  addition,  measurements 
made  with  a hot  wire  probe  behind  the  blade  row  indicated  that  substantial 
spanwise  flow  over  the  blades  existed  in  the  low  velocity  region  of  the  dis- 
tortion causing  the  flow  over  the  blades  to  be  three  dimensional.  Since  the 
analysis  holds  only  for  two-dimensional  flow,  only  qualitative  comparisons 
between  the  theoretical  and  experimental  results  can  be  made  although  in  each 
case  studied  the  analytical  upstream  axial  velocities  were  essentially  matched 
to  those  obtained  experimentally. 

Results  for  the  typical  case  corresponding  to  an  average  (3i  = 42.2°  are 
shown  in  Figs.  5 - 9.  Figure  5 shows  the  axial  velocity  of  distortion 
measured  approximately  1.2  chord  lengths  in  front  of  the  blade  row  leading  edge 
and  the  corresponding  distortion  predicted  analytically.  It  is  believed 
that  the  agreement  is  within  the  accuracy  of  the  experimental  results. 

A comparison  of  the  static  pressure  rise  across  the  blade  row  as  de- 
termined experimentally  and  analytically  is  presented  in  Fig.  6.  It  is  seen 
that  the  analytical  curve  exhibits  the  same  general  behavior  as  the  experimental 
data.  The  differences  in  the  magnitudes  of  the  two  results  are  thought  to  be 
due  to  errors  in  the  experimentally  determined  loss  curve  and  exit  flow 
angle  curve  used  in  the  theoretical  model  and  the  fact  that  the  flow  through 
the  blade  row  was  not  two  dimensional.  It  should  be  noted  that  changing  the 
magnitude  of  t,  the  time  constant  of  the  boundary  layer  time  delay,  in  eq.  (22) 
causes  the  analytical  pressure  rise  curve  to  shift  in  the  circumferential 
direction  without  any  appreciable  change  in  its  shape. 

Figures  7-9  are  "carpet  plots"  of  the  static  pressure,  resultant 
velocity,  and  total  pressure  throughout  the  entire  analytically  defined  flow 
field.  The  slight  "ripple"  in  the  axial  direction  of  the  total  velocity  and 
hence  total  pressure  is  apparently  due  to  numerical  instabilities  in  the 
numerical  solution.  It  should  be  noted  that  since  the  plot  routine  reduced 
the  axial  distance  by  a factor  of  two,  these  instabilities  are  not  as  severe 
as  they  appear  in  the  figures. 
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Appendix 

The  coordinate  system  used  in  the  calculations  covers  the  entire  domain 
of  interest  with  a rectangular  grid  of  cells,  each  of  dimensions  Ax  by  Ay. 

The  cells  are  numbered  by  indices  i and  j,  with  i counting  the  columns  in 
the  x-direction  and  j counting  the  rows  in  the  y-direction.  The  field- 
variable  values  describing  the  flow  field  are  directly  associated  with  these 
cells.  Their  points  of  definition,  relative  to  a cell,  are  shown  in  the 
figure  below: 


Field-variable  layout 

As  pointed  out  if  Ref.  6,  the  placement  of  the  field-variable  quantities 
relative  to  the  mesh  is  of  considerable  importance  to  the  matter  of  conser- 
vation. In  fact,  no  other  arrangement  other  than  the  one  shown  in  the 
figure  appears  to  be  workable. 

The  conservation  equations,  blade  row  coupling  equations  and  associ- 
ated boundary  conditions  in  finite  form  are  as  follows: 

I.  The  Continuity  Equation: 


Ax^Ui+*sj 


ui-%j^  + Ay^vij+*s  ~Vij-^ 


(A.  1 ) 
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II.  The  Momentum  Equations: 


ui+ij  = ff^i+lj  “ Dij)  + Hy^i+Jsj-1  + Ui+J5j)  (vi+lj-h 


+ vij->s)  ' (ui+y+l  + ui+y)  (vi+lj+>s  + vij+*s)] 


+ lkC(ui+%j  + Ui-Jsj)2  - (ui+*j  + ui+>sJ)21  + H(pij 


n.  ) + J^t,  + u.  , .)  + M(u.  , . •, 

pi+l j ^ ax?v  i+£j  -[-hj1  Ay^v  i+hj+1 


♦ ui+y-i* + n-2(^  * vy 


(A. 2) 


'lj+%"  3Ay^Di j+1  " Dij)  + 4^tvij-%  + vij+*5)2  " (vij+f+  Vij+*5)2] 


+ + ui-y)  (vi-lj+h  + vij+^  " (ui+hj+l 


+ ui+*sj^  ^vi+lj+3s  + vij+3s^  + Ay^pij  " pij+l^ 


+ aSP<v1+1  j+Jj  + VljW  + Ay^Vij+£  + VijV 

+ n_2(Mt.  + Mt»-|u 
*■  ^'axz  Ay 2 Jui+Jsj 


(A. 3) 


The  Pressure  Equations 


n = r Ax2Ay2 

pi j L2(AX2+Ay2 


^Rij  + ^2(Ax^+Ay2 


(pi+lj  + pi-l^ 


+l-2(Ax2+Ay2)^  ^pij+l  + pi  j-1  ^ 


(A. 4) 


where. 


"ij  - - 3ax^'  l+l j T ui-lj'  ‘ 3aP  v ij+1  i j-1 


* t§<ix2  4 ay2*  ' st]01j  + 44P[(,ij^  * vij+*5)3 

* (VU4  * W2  ' 2(vij«.  * Wij-%|J] 

+ 2fixay^ui+y+l  + ui+y^  ^ vi+1  j+*s  + vij+*s^ 

* (ui-y-l  * Vy*  + V1  j-’s’  ' (uf+>sJ-1 

+ Vy>  (Vij-b  4 W - <Vy+l 

4 Vy’  <V1-W  4 #W]  4 JSPt(ui**i 


4 ui+y>2  4 <V*j  4 ui-y>2  - 2(Vy  4 Vy,2] 


The  Blade  Row  Coupling  Equations: 


Pij  ■ P,+3j  4 4 <y  - Vy  - Vy> 


V c sec  c 


W~ vui+y+l  ui+%j  ui -W  i-*sj+l 

" V>sj-1  " Ui  +*sj  ' ui-*sj-l  ' Ui -*5J ) + ^Ui+bj 


4 ui-y,2[xj  4 (xj'1)  cot2  6u  4 C0t2  B2j] 


i+3j 


c sec  £ 
2At 


<ui+y  4 “i-y 


- U 


n+1 


i+y  " Ui-15j ) 


n+1 


+ n. 


ij 


+ 4Ay  ^ui+%j+l  + ui-hj+l  ~Ui+*ij-l  ~Ui -*sj-l 
- *<«<«,<  * * (*.  - 1)  cot2  B„  tcot2B2il 
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eij  ■ tan'’r(Vy  * Vy>  ' (2V  * * 'ij-ii11 

(A. 8) 

s2j  “ tan-’tCu^j  * ujtfj)  / (2V  + vj+3J+ls  * *„_*>] 

(A. 9) 

*?'  ■ “ *»*"  ■&  h * &vi  - xm> 

J 

(A. 10) 

Note  if  Xn+1  = X.: 

J J 

Xj  = xss.  + 2^Xj+l  " xj-l} 

J 

(A. 11) 

And  if  Xj+,  - ly 

XJ  ■ Xss  , 

J 

(A. 12) 

Inlet  Boundary  Conditions 

,i-U->5  “ Vi,j-H  • ’ v«+% 

(A. 13) 

D1-l j " °i j 

(A. 14) 

■c  - * w - 'C + O2 

* <Vy  + Vy>2  - <u"ty  * u"-y)2]  + pij 

(A. 15) 

u . ■» . = - u . .1  . + 2u . , • 
l-fj  i+’sj  i -*sJ 

(A. 16) 

Outlet  Boundary  Conditions 

Vi+lj+*5  = vij+*s  vi+lj-Js  = vij-*5 

(A. 17) 

Di+lj  ” °ij 

(A. 18) 

ui+-fj  ” ’ui-5sj  + 2ui+*sj 

(A. 19) 

VIII.  Blade  Row  Exit  Boundary  Conditions: 
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pi-l j = Ay^'Vy+l  ‘ Ui-*sj-l  + 2ui-y) 

+ 4Ay^ui-Jsj+l  + ui -*sj ^ ^ vi j+>s  + vi-lj+J^ 

- + W 1^+  + vi-ij-Ss)^ 

+ ftK-lj  " ui-y*  + i*lui+y  + Vy)2 

- (ui4j  + ui-^)2]  + pij  (A-28) 


Representation  of  Blade  Row. 


Observer  Reference  Frame 


Cascade  Reference  Frame 


Figure  2 Velocity  Diagrams  in  Stationary  and  Rotating  Reference  Frames. 


Upstream  Axial  Velocity  Distortion  Profil 
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